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The behavior of several nuclear properties with temperature is analyzed within the framework 
of the Finite Temperature Hartree-Fock-Bogoliubov (FTHFB) theory with the Gogny force and 
large configuration spaces. Thermal shape fluctuations in the quadrupole degree of freedom, around 
the mean field solution, are taken into account with the Landau prescription. As representative 
examples the nuclei ^®*Er, ^^^Dy and ^^^Hg are studied. Numerical results for the superfiuid to 
normal and deformed to spherical shape transitions are presented. We found a substantial effect of 
the fluctuations on the average value of several observables. In particular, we get a decrease in the 
critical temperature (Tc) for the shape transition as compared with the plain FTHFB prediction as 
well as a washing out of the shape transition signatures. The new values of Tc are closer to the ones 
found in Strutinsky calculations and with the Pairing Plus Quadrupole model Hamiltonian. 

PACS numbers: 21.60.Jz, 21.10.Ky, 21.10.Ma 



Introduction 

Since the advent of the new generation of 47r gamma 
ray detectors and the improved accuracy in the channel 
selection new possibilities have opened up in the study of 
nuclear structure. Besides this, the availability of faster 
computers has made possible to perform realistic theoret- 
ical investigations with large configuration spaces. The 
high excitation energy is specially interesting since new 
features may take place. For example, in the quasicontin- 
uum, the high level density gives rise to the unexpected 
phenomenon of the damping of the rotational motion. In 
the limit of high excitation energies (or temperature T) 
quantum effects become less relevant or may even disap- 
pear. Thus one expects that in a heated nucleus phys- 
ical effects like superfluidity or shape deformations are 
washed out when T increases. This expectation can be 
easily understood in terms of the shell model since, by 
increasing T, one promotes particles from levels below 
the Fermi surface to levels above it. In the case of pair- 
ing correlations, blocking levels amounts to destroying 
Cooper pairs. In the case of shape deformation, by de- 
populating the deformation driving levels (intruders) one 
gets on the average less deformation. Experimental in- 
formation about nuclear shape changes can be obtained 
by means of the Giant Dipole Resonance (GDR) built 
on excited states. Exclusive experiments studying the 
GDR strength at a given excitation energWor T) of the 
nucleus have been carried out in refs. 

Hlilill The 

understanding of these phenomena is relevant because it 
affects important features like the fission barriers and the 
stability of the nucleus itself. For a recent review on hot 
nuclei see ref. (5|. 

The shape transitions have been object of many stud- 
ies, most of them with schematic models, separable 
forces, and small configuration spaces 0, 0, BM tl0| - 
The theoretical approaches used in the calculations are 



based on the mean field approximation, mainly the Finite 
Temperature Hartree-Fock-Bogoliubov theory (FTHFB) . 
The mean field approximations predict sharp shape tran- 
sitions, whereas for finite systems, however, one expects 
washed out transitions instead. The fact that the pre- 
dicted critical temperatures are rather high (around 2- 
3 MeV) indicates that not only the most probable de- 
formation is relevant but that there is a finite (in some 
cases very large) probability for the system to have other 
shapes which should be taken into account. Calcula- 
tions beyond mean field including thermal fiuctuations 
have confirmed the expectation of washed out transitions 

iniiia. 

Theoretical studies with effective forces and large con- 
figuration spaces have been performed at the FTHFB 
level with density dependent forces, Skyrme 0, 0, 
and recently with the Gogny force. Additional cal- 
culations have been done in the relativistic mean field 
(RMF) approximation [13, 0|- Calculations including 
thermal fluctuations in conjunction with large configu- 
rations spaces and effective forces have been performed 
only very recently fvA . 

From these studies a discrepancy has emerged since, 
while the mean field approaches (FTHFB) with effective 
forces (like Skyrme, the Gogny force or the relativistic 
approaches), provide the view of a sharp shape transi- 
tion at a relatively high critical temperature (Tc ~ 2.7 
MeV for ^^'^Er), schematic models (like the Pairing plus 
Quadrupole) and Strutinsky calculations provide also a 
sharp transition though at a much lower critical temper- 
ature {Tc « 1.7 MeV for ^^''Er). Furthermore, the discor- 
dant point of a "sharp" transition for a small system, like 
the nucleus, predicted by both approaches requires fur- 
ther investigation. Earlier calculations with the Pairing 
plus Quadrupole Hamiltonian ^ll, 20J have pointed out 
the relevance of including fluctuations in mean field ap- 
proaches at finite temperature as a step forward to clarify 
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some aspects of these problems. It is the aim of this pa- 
per to investigate the problems just mentioned as well as 
other related high excitation energy topics, level densi- 
ties, etc., within a beyond mean field theory. Towards 
this end, the FTHFB calculations of ref. JS} with the 
Gogny force and large configuration spaces will be gener- 
alized to include fluctuations in the quadrupole moment 
degree of freedom. 

The finite range density dependent Gogny force has the 
advantage of providing the particle-hole (Hartree-Fock) 
and the particle-particle (pairing) matrix elements from 
the same interaction, at variance with relativistic theo- 
ries and most Skyrme calculations. In the fitting of the 
DIS |23| parametrization no excited states or spin de- 
pendent data was used, however, it has produced good 
results in the description of nuclear properties not only 
at zero Ip^. l23l | but also at large spin |3| and, more 
recently [l6j, in calculations at high excitation energy. 
Since our purpose is to study the behavior of shell ef- 
fects and fluctuations with temperature, we have selected 
both theoretically and experimentally well known nu- 
clei that display a variety of shapes in the ground state: 
strongly deformed (^^^Er), oblate (^^^Hg ) and rather 
soft i^^^By). 



Theory 

To study the behavior of nuclei with increasing tem- 
perature we use the DIS parametrization of the fi- 
nite range density dependent Gogny force |22i, in the 
FTHFB framework 0, The Gogny force, at vari- 

ance with most of the Skyrme parametrizations and the 
relativistic models, allows full selfconsistent calculations 
since it provides the particle-hole and pairing fields from 
the same force. 

At finite temperature, as at temperature zero, the ba- 
sic approximation is the mean field theory. Its most so- 
phisticated version, the FTHFB, has been developed in 
refs. 0, S ll^ . For convenience we will give here a short 
outline. 

For a system at constant temperature T and with 
chemical potential fj,, the equilibrium state can be ob- 
tained from the variational principle over the grand 
canonical potential 



operator given by 



E ~TS ~ fiN. 



(1) 



The energy, E, entropy, S, and particle number, N, are 
thermal averages defined by 

E = (h^^^ = Tr{DH), 

-fclnZ)^^ = --kTr{D\nD), (2) 
N = (^n"^^ = Tr{DN), 

with Z the grand partition function and D the density 



Z = rr[cxp(-/3(H-M^))], 
D = Z-^cxp{-/3{H - fiN)), 



(3) 



with P = 1/kT. 

In the FTHFB approach the density operator is ap- 
proximated by 



Do 



(4) 



where Ti. is the most general Hermitian single particle 
operator, to be determined by the variational principle 
and Zo is the partition function. It can be shown 5] that 
TC is given by 



n 



h A 

-A* -h* 



(5) 



with A the pair potential and h the HF hamiltonian. h 
is given in terms of the kinetic energy t, the HF field, F, 
and the chemical potential, /i. 



h — t + r — /i. 



(6) 



The density matrix, p, and the pairing tensor, k, are 
given by 



and 



p = UfU+ + V*il- f)V\ 
K ^ UfV+ + V*{1- f)U\ 



1 + e/^s. 



(7) 



(8) 



The matrices (C/, V) provide the relation between the 
quasiparticle and the single particle basis : 



fe 



(9) 



They are determined, together with the quasiparticle en- 
ergies, Ei, by the FTHFB equation 



h A 

-A* -h* 



Uk 
Vk 



Uk 
Vk 



Ek 



(10) 



The solution of this equations provides us with the con- 
figuration that minimizes the grand canonical potential. 
With U, V and / known one can determine the density 
operator Dq and calculate any expectation value. For 
density dependent forces like Skyrme or Gogny, the for- 
malism remains unchanged except in the evaluation of 
the one body Hamiltonian h. Due to the dependence on 
the density of the interaction, h gets [T9l| an extra term, 
dr, which is usually referred to as the "rearrangement 
potential" and is given by 



dVr. 



dH 



(11) 
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The FTHFB solution gives us the most probable 
shapes, quadrupole, hexadecupole, etc, as well as the 
most probable gap parameters and so on. At finite tem- 
peratures, however, we have statistical (or thermal) fluc- 
tuations around this solution. In principle one could con- 
sider fluctuations in the most relevant degrees of freedom. 
For nuclei, at high excitation energy, the most impor- 
tant one is the quadrupole deformation, and we therefore 
shall consider only the fluctuations in the quadrupole mo- 
ment {Q20) in this paper. To generate the solutions with 
different shapes we solve the grand canonical potential, 
Eq. with an additional constraint on the quadrupole 
moment, i.e., we minimize fl — E —T S — — XQ^^q. The 
Lagrange multiplier Aqjo is adjusted in such a way that 
the thermal expectation value {Q20) = Tr(DQ2o), has 
the required value q. According to Landau |27J the prob- 
ability P{q) to obtain a certain value q of the deformation 
is characterized by the free energy F{q) = E{q) — TS{q) 
of the system with deformation q 



P{q) oc e-^(«)/^. 



(12) 



Using classical statistics, therefore, for the ensemble av- 
erage of an observable O one obtains the expression 



O 



JO{q)eM~F{q)mdq 
JeM~Fiq)/T)dq ' 



(13) 



where 0{q) is the thermal expectation value of the op- 
erator O calculated for the system with the deformation 
(?, and dq is the volume element in deformation space. 
In our case the set q corresponds to the quadrupole de- 
formation q2o, thus dq = dq2o, with metric equal to one. 
The limits in the thermal average integrals, see Eq. p3|l . 
are chosen to span the full (32 region in which the prob- 
ability of having one of these values, given by Eq. (|12|l . 
is not negligible. This covers both prolate and oblate 
regions. 

High temperature calculations require large configu- 
ration spaces. In order to maintain the computational 
burden within reasonable limits we restrict ourselves to 
axial symmetry. We are aware that for soft nuclei and/or 
high temperature, triaxiality may play an important role. 
In the calculations we use an axially deformed harmonic 
oscillator (HO) basis with a size defined by the condition 



2hpnp + h^n^ < Nq, 



(14) 



where Up and Uz are the HO axial quantum numbers, 
bp = q^/'^ and bz = q~'^^^ with q = Rz/Rp, the nuclear 
axis ratio. In our case we have used q = 1.5 and A'o = 15 
which allows for deformations big enough to reach the 
fission barrier and provides room enough for the tem- 
perature induced excitations. However, as an additional 
check, we have also used A^o = 17 for some selected calcu- 
lations. Reflection asymmetry is allowed in the calcula- 
tions, i.e. the nuclei may develop octupole deformations. 

In order to compare our Gogny force results with the 
ones, more conventional and popular, of the schematic 



Pairing plus Quadrupole model (PPQ), we have also per- 
formed calculations with this force. The configuration 
space (the spherical oscillator shells N — 5,6 for neu- 
trons and A = 4, 5 for protons) and the force parameters 
used are the one of Baranger-Kumar _28]. The calcula- 
tions have been performed in exactly the same way as in 
ref. El 



Results 

We have performed FTHFB calculations with the DIS 
parameter set of the Gogny force in several nuclei to 
study the evolution of shell effects with temperature. Nu- 
clei with different ground state deformations have been 
selected to illustrate their different behavior. As an ex- 
ample of a nucleus with a strongly prolate deformed 
ground state we used the thoroughly studied ^^^Er. The 
soft ^^^Dyse is a transitional nucleus between the clearly 
spherical Dysprosium isotopes with A^ < 84 and the well 
deformed ones with A^ > 88. This nucleus was selected 
for its rich shell structure and shape coexistence. The 
heavier ^^^Hg has been chosen due to its oblate ground 
state shape. The maximum temperature studied has 
been kept below 3 MeV, such that continuum contribu- 
tions can be safely disregarded |2^ • 

The thermal fluctuations are represented through av- 
erages calculated according to Eq. (fT^ . The deviations 
around the mean values can be studied by the standard 
deviation value 



cr(O) = W 02 - O 



(15) 



This quantity is presented in some cases for further clar- 
ification of the results obtained. 

Before entering in the discussion of the shape and pair- 
ing phase transitions we will start by presenting first a 
general view. In Fig.^we present the free energy, F{/32), 
and the quantity P{P2) oc exp(— i^(/?2)/r) versus the 
quadrupole deformation, (32, at different temperatures. 
P{P2) provides the weight of a given shape (32 in the eval- 
uation of thermal average values. The results for ^^''Er 
are displayed in the left column, for ^^^Dy in the middle 
one and for ^^^Hg in the right one. For the rare earth 
nuclei -'^^''Er and ^^^Dy, where calculations within the two 
shells configuration space mentioned above are feasible, 
we also present results with the PPQ model. Results 
with the Gogny force are displayed by continuous lines 
and those with the PPQ force by dashed ones (thick lines 
represent F{(32) and thin ones P{(32))- The well depths 
are measured from the point with (32—0 and P{P2) has 
been normalized is such a way that the most probable 
deformation takes the value of unity. 

Let us discuss first the low temperature calculations 
(T = 0.3 MeV ) where we can observe the intrinsic shapes 
of the ground states. In the Gogny calculations for ^^"'Er 
there is a deep prolate minimum at (32 ~ 0.3 and about 
4.5 MeV higher an oblate one. With the PPQ model 
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the same gross features are observed, though the min- 
ima are not so deep. The probabihty distribution P(/32), 
however, is similar in both calculations. For the nucleus 
i52j3y^ in the Gogny case, the prolate minimum is at 
/32 ~ 0.15 and a bit higher in energy the oblate one. 
The PPQ model, for this nucleus, provides a broad mini- 
mum around the spherical shape. P(/32), in contrast with 
^^''Er, looks quite different in the Gogny case than in the 
PPQ one. In both nuclei the free energy surfaces are 
broader with the Gogny force than with the PPQ one. 
For ^^^Hg, we find the minimum at an oblate deforma- 
tion of (32 ~ —0.15 and about 1.7 MeV higher a small 
prolate minimum. 

At higher temperatures the expected disappearance of 
shell effects becomes clear, in particular the vanishing of 
the barriers when several minima are available and the 
development of only one spherical minimum. Further 
finite temperature effects like the widening of the free 
energy curve and the more important role of fluctuations 
with increasing temperature appear in the Gogny calcula- 
tions but not in the PPQ ones. In the PPQ case the free 
energy surfaces, with increasing temperatures, become 
flatter but not broader. They even become narrower! 
This unphysical effect has to do, obviously, with the size 
of the configuration space (two shells). As we can see 
already at (T = 0.3 MeV ), at large deformations F{/32) 
increases very steeply because there are no orbitals with 
high-j (deformation driving ) coming down. The mech- 
anism to soften the free energy surface at high temper- 
ature by enhancing the probability to occupy high-lying 
orbitals (among them the high-j ones) works only with 
large configuration spaces. The anomalous behavior of 
P{p2) in ^'■^^Dy, in the PPQ approach, at T = 0.6 MeV 
as compared with T = 0.3 MeV is due to the fact that at 
T = 0.6 the neutron pairing gap vanishes and the ground 
state becomes slightly prolate. 

Figures 121 El and El show the detailed calculations for all 
three nuclei. These figures include both, the results at the 
FTHFB level and with shape fluctuations calculated as 
described above. Dashed lines and open symbols indicate 
the FTHFB results. Solid lines and filled symbols are 
used for averaged, fiuctuations including, calculations. 
Figs. Eland El show only Gogny results. 



A. The nucleus "''Sr 

In Fig. 121 we display the results of the calculations for 
the nucleus ^^"'Er with the Gogny force and with the 
PPQ model Hamiltonian. In panel (a) we show the self- 
consistent FTHFB (i.e., calculated with the solution of 
Eg. I10|l and the averaged results (i.e., calculated accord- 
ing to Eq. I13|l for the (32 , Pi and (3q deformation param- 
eters as a function of the temperature with the Gogny 
interaction. Let us first discuss the deformation param- 
eter ^2- For temperatures < T < 1.0 MeV, both pre- 
dictions behave similarly, as one would expect for a nu- 
cleus with a well pronounced minimum. For tempera- 



tures 1.0 < T < 2.0 MeV, the FTHFB /?2-values de- 
crease rather smoothly while the averaged ones undergo 
a strong reduction. For T > 2.0 MeV the selfconsistent 
values decrease very steeply and collapse, finally, to zero 
deformation at T = 2.7 MeV. The averaged values, on 
the contrary, change tendency decreasing very smoothly 
in such a way that an almost constant value of /?2 is 
eventually obtained. The behavior of (3^ and (3q is sim- 
ilar to the one of /32 though not that spectacular. The 
same plot for P2 but with the PPQ interaction is repre- 
sented in panel (b). Quantitatively the main differences 
with the Gogny results are the faster collapse of the self- 
consistent value, at T « 1.8 MeV, and the reduction of 
the temperature interval where the averaged values are 
smaller than the selfconsistent ones. Looking at the prob- 
ability distribution in Fig. ^ one can easily understand 
these differences. The temperature value at which the 
mean field (mf) deformation parameter collapses, which 
we will denote T™^, has often been used in earlier mean 
field studies to signal a shape phase transition. It is obvi- 
ous from panels (a) and (b) that in theories beyond mean 
field things look quite different and that definition of the 
critical temperature must be carefully considered. The 
big difference in as predicted by effective forces, like 
the Gogny force, and the PPQ is also known from cal- 
culations with Skyrme forces I13I, IT^ and the relativistic 
mean field approximation |l7l l3l|. 

The standard deviation in the deformation parameter 
f32, f (/32), calculated according to Eq.^l is presented in 
panel (c). One can distinguish three well defined zones : 
in the first one at low temperature, when pairing is still 
strong, the deformation is kept almost constant and fluc- 
tuations raise slowly. The Gogny and PPQ calculations 
predict about the same equilibrium shape in this zone. 
For temperatures higher than the one corresponding to 
the pairing collapse, cr(/?2) increases rapidly up to a maxi- 
mum value, remaining more or less at this value at higher 
temperatures. This step behavior is characteristic of a 
shape phase transition region. In fact, one could deflne 
the shape transition temperature as the one at which 
cr(/32) has a maximum. With this criterion one obtains 
T ^ 1.4 MeV for the PPQ result and T ^ 1.7 MeV in the 
Gogny case. Note the large difference in a {(32) between 
the Gogny and PPQ results at high temperature. In this 
comparison, however, one must keep in mind that the 
PPQ model hamiltonian is restricted to a configuration 
space of two oscillator shells which strongly constraints 
the ability to produce fluctuations. This is clearly seen 
in Fig. n : the PPQ results rapidly develop a narrow 
parabolic shape with increasing temperature. This lack 
of fluctuations was already identified as partially respon- 
sible for the low multiplicity seen in the collective E2 qua- 
sicontinuum spectra in gamma decay calculations when 
compared to experiment jHM 1^ ■ 

An additional confirmation of the importance of the 
fluctuations in calculations with the Gogny interaction as 
compared to the PPQ case is provided by the different 
behavior of the averaged deformation parameters with 
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respect to the mean field within the same model, see 
panels a) and b). The value of the average /?2 parameter 
at high temperatures illustrates the deviation of the free 
energy surface from a parabolic behavior. A value close 
to zero is expected for a parabola, e.g., in the PPQ case, 
while a larger one, as in the Gogny case, indicates the 
softness of the prolate side as compared with the oblate 
one. One should nevertheless keep in mind that only 
axially symmetric deformed shapes are allowed in the 
calculations. 

In panel d) of Fig. [3 the proton and neutron pair- 
ing energies are displayed for the Gogny force. Up to 
T « 0.2 MeV the pairing energies are rather constant 
but for higher T values they decrease in absolute value 
very fast up to T = 0.7 MeV where they vanish. Ther- 
mal shape fluctuations, as expected in the low tempera- 
ture regime, have little effect on the pairing correlations. 
Pairing fluctuations which would be more relevant [2(1 | 
are not considered in this work. 

It is interesting to take a look at the internal excitation 
energy, E* , and also analyze it through the behavior of 
its derivative, the specific heat Cv{T) — dE* /dT, since 
the appearance of peaks in this quantity is customarily 
interpreted as a signature in the search for phase transi- 
tions. 

The evolution of E* with temperature for both the 
FTHFB, E* , and the average, E*, calculations, is pre- 
sented in panels e) (Gogny force) and f) (PPQ force) of 
Fig. 12 With the Gogny force and in the low temperature 
regime, we can see the pairing collapse which is visible 
as a change in the slope of E* . At higher temperatures a 
fairly quadratic behavior is observed in the excitation en- 
ergy, which is slightly modified when the transition to the 
spherical phase takes place at high temperature around 
T ~ 2.7 MeV. There, a weaker change in the slope, hardly 
seen in the scale of the plot, is found. The change is more 
abrupt in E* than in E* , again as expected in the pic- 
ture of a thermally faded transition. The same facts are 
observed in the PPQ plot, where one can additionally ob- 
serve that at high temperature the energy behaves more 
linear than quadratic as a function of T. 

The different behavior of E* and E* is also interesting : 
a) At temperatures below 0.8 MeV, both energy values 
coincide, b) between 0.8 MeV and the corresponding 
(around 2.7 MeV for Gogny and 1.8 MeV for PPQ), E* is 
always below E* , c) at the critical temperature both en- 
ergies do coincide and d) at higher temperatures they are 
rather similar. This behavior has a simple explanation if 
one considers the entropy as a function of the deforma- 
tion at fixed temperature, see Fig.|3| and the fact that in 
general F < F ,36] , with F the selfconsistent FTHFB free 
energy at the given T. At low temperatures (< 0.8 MeV) 
and at temperatures above T™^, the entropy is almost 
shape independent, i.e., S S (with S the selfconsistent 
FTHFB entropy at the given T) , and the free energy be- 
haves like a parabola, F « F, consequently E* w E*. 
At temperatures 0.8 < T < T^"^, F < F_and _S < 'S, 
see Figs, n] and 13 and consequently, since F = E — TS, 



E* < E*. 

The change in slope in E* as compared with E* indi- 
cates that the corresponding specific heats, dE* /dT and 
dE* /dT, will be rather different. This can be seen in 
Fig. 21 where we observe two peaks in the selfconsistent 
results, both in the PPQ and in the Gogny calculations 
|37l |. In between we find the typical linear behavior for 
a Fermi gas (Cy — 2aT). The low temperature peak is 
associated with the superfluid to normal transition and 
the high temperature one with the deformed to spheri- 
cal shape transition. The low temperature peak remains 
nearly unaffected by the inclusion of fluctuations. At the 
mean field level this transition takes place at a tempera- 
ture low enough such that shape fluctuations are irrele- 
vant. 

Comparing the FTHFB results with the PPQ force 
with those obtained with the Gogny interaction for the 
shape transition, we find again the different temperature 
predictions. Using as critical temperature, Tc, the tem- 
perature where Cy changes curvature, the same values 
are obtained as when the (32 — ^ rule was used (1.8 MeV 
for PPQ vs. 2.7 MeV for Gogny). However, when fluctu- 
ations are taken into account, the difference in the pre- 
dicted Tc by Gogny and PPQ gets smaller and the sharp- 
ness of the peaks reduced, indicating a less abrupt transi- 
tion as is expected in a mesoscopic system. This is more 
evident in the Gogny results, where the peak becomes a 
broad bump, providing another clue of the greater im- 
portance of fluctuations in the Gogny case. By contrast, 
the PPQ peak, although broader than in the mean field 
case, is still sharp. Furthermore the PPQ specific heat 
levels off, showing how the limited configuration space 
available is a clear disadvantage of this model. If we now 
look for the changes in curvature we find Tc — lA MeV 
in the PPQ and Tc = 1.7 MeV with the Gogny interac- 
tion. It is interesting to notice that these values agree 
very well with the ones obtained looking at <j{P2)- 



B. The nucleus ^^^Dy 

As we have seen in Fig. ^ ^^"^Dy displays a potential 
energy surface with energetically close prolate and oblate 
minima. It could illustrate a nucleus with shape coexis- 
tence, that means, already at temperatures near to zero 
there is a finite probability of populating more than one 
minimum. For this nucleus we will not perform a discus- 
sion as exhaustive as for ^^'^Er but we will consider the 
most relevant facts. 

In Fig. [SI the results of the calculations for the Gogny 
force are displayed. In panel a) the mean field and the 
averaged values of the pairing energies are plotted. The 
pairing energy of the neutron (proton) system collapses 
at T « 0.5 (T « 1.0) MeV. The averaged values, as 
expected, almost coincide with the mean field ones. In 
panel b) the /3-deformation parameters are shown. As we 
can observe in the behavior of the [32 parameter the effect 
of the shape fiuctuations in this case is already noticeable 
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at very small temperatures. This is due to the fact that, 
in the ground state ^^^Dy is a much less deformed nu- 
cleus than -'^^'^Er and that the energy difference between 
the oblate and prolate minimum is small, amounting to 
only ^ 0.55 MeV for T < 0.5 MeV. Hence the averaging 
formula assigns finite weights to the oblate side already 
at small T's causing the observed steep decrease in the 
average P2 value. Although the FTHFB, searching for 
the strict minimum, provides a deformed ground state 
for ^^^Dy, we see how even at the lowest temperatures 
the average deformation is very small, in agreement with 
the experimental data. At the high temperature limit 
we observe that above 1.4 MeV, the average deforma- 
tion stays rather constant, or slightly increases, up to 
0.04. This anomalous behavior is due to the fact that 
the superdeformation driving orbitals are being occupied 
at this temperature range. The /34 and Pe deformation 
parameters follow closely the P2 behavior. In particular, 
in the FTHFB description they become zero at the same 
temperature as /?2. 

In panel c) the excitation energy is depicted. In the 
FTHFB approach the changes in slope at temperatures 
of 0.5, 1.0 and 1.4 MeV are due to the neutron and pro- 
ton pairing collapse and to the shape transition. The be- 
havior of E* can be understood in terms of the entropy 
plots as explained for the ^^''Er nucleus. The changes 
in slope in the energy plots are magnified in the spe- 
cific heat versus temperature plot depicted in panel d). 
In the mean field approach we find a broad composite 
peak, corresponding to the proton and neutron pairing 
collapse, showing substructures around T « 0.5 and 1.0 
MeV. Furthermore, one sees a second peak at 1.4 MeV 
corresponding to the FTHFB shape transition. If shape 
fluctuations are included in the calculations we obtain 
only one broad peak. The small peak at 1.4 MeV, how- 
ever, is not there anymore indicating that it has been 
shifted to the pairing transition bump or simply washed 
out. In fact, the almost identical broad energy mean field 
peak and the single one obtained with fluctuations could 
be seen as a clue that the shape transition in this soft nu- 
cleus is inexistent, since there is no higher temperature 
peak but only a small modification of the "pairing" one. 
To check this hypothesis we have also performed calcula- 
tions with the particle-particle channel of the Gogny force 
set to zero, i.e. plain FTHF. In this way we obtain Cy 
curves without the pairing transition peaks. The results 
are plotted, superimposed to the standard FTHFB calcu- 
lations, using square symbols in panel d) of Fig.|Sl Again, 
open squares are for the FTHF calculations and full ones 
for averaged ones. As it was expected, the FTHF curves 
for Cv show no peak for the pairing transition and the 
only peak present is the one corresponding to the shape 
transition at 1.4 MeV, which coincides with the one ob- 
tained with the full FTHFB calculations, since at these 
temperatures pairing is already zero. The thermal aver- 
aging results show a broad shoulder approximately in the 
same temperature region in which the pairing transition 
was located. This rather soft bump is a clear indication 



of the above mentioned situation, i.e., at small temper- 
atures no clearly predominant minimum exists and at 
high temperatures the nucleus does not become exactly 
spherical. 



C. The nucleus ^^^Hg 

As we have seen in Fig. ^ at very small temperatures 
this nucleus presents an oblate deformed ground state 
and about 1.7 MeV above a prolate minimum. The re- 
sults for the nucleus ^^^Hg are shown in Fig. El The 
pairing energies are displayed in panel a) . The proximity 
of the Z = 82 shell closure causes the vanishing of the 
proton pairing energies for all temperatures. The neu- 
tron system, on the contrary, has a large pairing energy 
at T = MeV, which vanishes at T = 0.8 MeV. As be- 
fore, the shape fluctuations have almost no effect on the 
pairing energies. 

In panel b) the behavior of the deformation parameters 
with increasing temperature is plotted. For ^^^Hg, in the 
mean field approximation, we obtain an oblate ground 
state deformation of /32 = —0.135 which gets more oblate 
for increasing temperatures as the pairing energies go to 
zero. For larger T values the deformation decreases and 
around T « 1.4 MeV the nucleus becomes spherical. As 
before, the effect of the shape fiuctuations is mainly char- 
acterized by the the prolate-oblate ground state energy 
difference, see Fig. ^ which in this case amounts to 1.7 
MeV at zero temperature and only above T = 0.5 MeV 
start to diminish. Around this temperature the defor- 
mation gets smaller and around 1.4 MeV the average 
deformation is zero. Interestingly, the average P2 and (3^ 
deformation parameters become positive in the limit of 
high temperatures due to the fact that at these temper- 
atures the F{T) curves are softer in the prolate than in 
the oblate side. 

In panel c) the excitation energies in both approxima- 
tions, in the mean field and with shape fiuctuations, are 
plotted versus the temperature. In the mean field calcu- 
lations we find slope changes at T = 0.8 MeV and T = 1.4 
MeV associated with the pairing collapse and the oblate- 
spherical shape transition. The inclusion of shape fiuc- 
tuations affects mainly the region between 0.8 <T< 1.4 
MeV. The general behavior of both curves can be easily 
understood in the same terms as for -'^^''Er. In panel d) 
the corresponding specific heats are represented. In the 
mean field approximation the expected peaks are clearly 
visible. The inclusion of the shape fluctuations produces 
a single, broader bump extending above the critical tem- 
perature for the pairing collapse. As in the ^^^Dy case 
this peak might be a superposition of the pairing and 
the shape transition peaks. To isolate the shape transi- 
tion peak we have performed again calculations with the 
particle-particle channel of the Gogny force set to zero. 
The results of the calculations, without (empty squares) 
and with (filled squares) shape fluctuations, are repre- 
sented in the same panel. In both approximations the 
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results above T = 0.9 MeV are obviously the same as 
before. Below this temperature and in the mean field 
approximation, as expected, the pairing peak is gone. 
However, in the calculations with the shape fluctuations 
we find a broad peak extending from T — 0.4 MeV up to 
T = 1.0 MeV with a change in curvature around T = 0.9 
MeV. Looking at the standard deviation o'(/?2) of this 
nucleus (not shown here) we find a maximum at T = 0.9 
MeV, an additional indication of the shape transition. 



D. Level densities and nuclear radii 

Level densities, p{E*), can be microscopically evalu- 
ated in the saddle point approximation, see for example 
eq. (2B-14) of ref. >3l|. In Fig. □ the total level densities 
for the three nuclei under study are displayed against the 
excitation energy in the mean field approximation and on 
average, i.e., with the inclusion of shape fluctuations. In 
both cases we observe the overall expected exponential 
dependence and the well known abnormal behavior at 
very small excitation energies. 

For ^'^''Er and up to 10 MeV excitation energy (T « 0.8 
MeV) we find a good agreement between both predic- 
tions. Then, up to 70 MeV (T « 2.3 MeV), we observe 
an increase in the level density in the average descrip- 
tion as compared with the mean field one. In particular, 
around 20 to 30 MeV excitation energy, the average pre- 
scription provides almost two orders of magnitude larger 
densities than the mean field one. This behavior can be 
easily understood looking at Figs. ^ and O and taking 
into account that the level density is proportional to the 
exponential of the entropy. At low (T < 0.7 MeV) and 
high excitation energies (T > 2.5 MeV) the entropy is 
rather shape independent, that means, the average value 
of the level density is very close to the one in the selfcon- 
sistent minimum. Consider now T = 1.4 MeV. In this 
case the selfconsistent minimum is prolate {(32 ~ 0.3), 
and the entropy at this shape and T is smaller than for 
all the other shapes at this temperature. That means, 
since p{(32) (X e^^'^^\ the average level density will always 
be larger than the selfconsistent one. Similar arguments 
apply to understand the behavior of the level densities 
of ^^^Dy and ^^^Hg. The fact that in these nuclei we do 
not find a larger difference between both descriptions is 
obviously due to the smoother behavior of the entropy 
with the deformation at the relevant temperatures. 

In Fig.|Slthe root mean squared (rms) radii of the three 
nuclei are plotted versus the temperature. In general we 
find that the rms radii are rather constant up to a given 
temperature, 2 MeV for ^^^Er and 1 MeV for ^^^Dy and 
^^^Hg, and that in this temperature range the average 
values are rather similar to the FTHFB ones. From this 
temperature on the average values are larger than the 
FTHFB ones due to the fact that, at these temperatures, 
the probability for a given shape peaks at the spherical 
shape and that for a given volume the spherical shape 
corresponds to the one with the smallest rms radii. That 



means, fluctuations around the spherical minimum pro- 
vide always larger rms radii. We also observe, at the 
highest temperatures, the expected increase of the rms 
radii. 



Discussion and conclusions 

We have seen in the previous section that, in calculat- 
ing average properties, the behavior of the entropy with 
the deformation parameter (32 plays a major role. This 
behavior is by itself, indeed, quite interesting. The gen- 
eral behavior, see Fig. O is the following : In the high T 
limit where the temperature effects dominate, we find, as 
expected, small shape dependence. At very low T, where 
the temperature effects are very small, we observe that 
to increase the entropy by 5 units we have to increase T 
by 0.5 MeV and that the entropy is rather independent 
of the shape of the nucleus. Of course, in this region 
where pairing correlations are present it is difficult to 
make more precise statements. However, at moderate 
temperatures, which are however high enough to allow 
for significant quasiparticle occupation numbers but not 
too high in order that shell effects are still present, one 
can find a large dependence of the entropy on the nu- 
clear shape. In this region spherical shapes, as expected, 
have larger entropy than axially deformed ones. Since the 
maxima of the entropy are associated with the minima 
of the grand potential a correspondence between Fig. Q] 
and Fig. 131 does not necessarily exist. 

One of the main outcomes of our research is the find- 
ing that shape fluctuations have a large effect on the de- 
scription of shape transitions. In fact, the mean field 
approach (FTHFB) with effective forces (like Skyrme, 
the Gogny force or the relativistic approaches), provides 
the view of a sharp shape transition at a relatively high 
critical temperature (T^ « 2.7 MeV for ^^'^Er). On the 
other hand, Strutinsky calculations or schematic models 
(like the Pairing plus Quadrupole) provide also a sharp 
transition though at a much lower critical temperature 
(Tc « 1.7 MeV for i^-^Er). It has been argued [II that 
the different predictions for the critical temperature are 
due to the small effective mass obtained in the mean field 
approach with effective forces {m* /m ^ 0.7, with m the 
nucleon bare mass ) as compared to the Strutinsky or 
the PPQ model (m*/m « 1.0). This argument is obvi- 
ously restricted to the mean field approach. In theories 
beyond mean field it does not apply anymore because 
with increasing correlations the effective mass eventually 
becomes the bare mass. In calculations at finite temper- 
atures two kinds of correlations have to be considered, on 
one hand the quantum ones and on the other the statis- 
tical or classical ones. Their relevance depends obviously 
on the excitation energy (or temperature), at low T's the 
former are very important and at high T's, the latter 
ones. At the temperatures where the shape transition is 
predicted to take place in the mean field approach, the 
probability of having a shape different to the selfconsis- 
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tent one is very large. Therefore, it is obvious that, first, 
shape fluctuations must be included and second that the 
characterization of the shape transition must be consid- 
ered more carefully. In the mean field approach a crite- 
rion for shape transition is just to look at the temperature 
at which the nucleus becomes spherical or alternatively 
to look for a peak at the specific heat as a function of the 
temperature. In theories beyond mean field, usually the 
second one is used because the average deformation can 
become very small but not zero. As we have seen, the 
inclusion of shape fluctuations provides a specific heat 
rather different from the mean field one, because a) it 
is not a sharp peak what we obtain but a rather broad 
bump (this is consistent with the fact that it is a very 
small system) and b) it appears at temperatures much 
lower than the ones predicted by the mean field approxi- 
mation in agreement with the Strutinsky calculations. It 
is also interesting to notice that the predictions based on 
the specific heat analysis coincide with the ones of the 
standard deviation (t(/32). 

Of course one could ask about the effect of considering 
quantum correlations in our predictions. Canosa, Rossig- 
noli and Ring |35j | have shown in model calculations based 
on the static path plus random phase approximation that 
at finite temperature quantum effects are observable de- 
pendent. In particular, they find that the specific heat re- 
mains practically unaffected when quantum correlations 
are taken into account. One could conclude therefore 
that the prediction of the shape transition by the inclu- 
sion of thermal shape fluctuations with effective forces is 



reliable. 

In conclusion, we have performed extensive calcula- 
tions with the Gogny force and a large configuration 
space for three representative nuclei in the FTHFB 
framework. We have further studied the effect of thermal 
shape fluctuations and found that they strongly affect, 
among others, the traditional shape transition "view" of 
the FTHFB approach. They do it in two aspects, first 
the critical temperature for the transition is very much 
lowered and second, the specific heat peaks are not sharp 
but rather broad. Besides this, the peaks showing up in 
FTHFB calculations of the specific heat in soft nuclei, like 
^^^Dy, are (almost) completely washed out when thermal 
shape fluctuations are taken into account indicating the 
absence of any shape transition. However, in strongly 
deformed nuclei, like ^^^Er, the shape fluctuations re- 
confirm the presence of a shape transition though of a 
different character. We also find a strong enhancement 
in the level density in the presence of a shape transition. 
The superfluid to normal phase transition is not affected 
by the inclusion of shape fluctuations. 
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FIG. 1: Free energy curves for ^''*Er, ^^^Dy and ^^^Hg calculated with the Gogny (thick solid lines) and PPQ (thick dashed 
lines) interactions at several temperatures as a function of the quadrupole deformation parameter /32- The probabilities P{P2) 
for a given shape /32, with the Gogny (thin solid lines) and the PPQ (thin dashed lines) interaction. The scale for P{fi2) is 
given in the inset. 
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FIG. 2: Results for ^®*Er for several observables as a function of the temperature. In all panels except in (d), solid lines and 
filled symbols are average values calculated according to Eg. 1131 dashed lines and open symbols are for selfconsistent (FTHFB) 
results. Circles are for Gogny and squares for PPQ. Shown are the (32 deformation parameter with Gogny (panel a) and PPQ 
(panel b) forces, the standard deviation (t(/32) with Gogny and PPQ forces (panel c), the pairing energies with the Gogny force 
(panel d) and the excitation energy with the Gogny (panel e) and the PPQ force (panel f). 
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FIG. 3: The entropy at fixed temperatures as a function of the deformation parameter ^2- 
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FIG. 4: The specific heat for tfie nucleus ^^''Er with the Gogny (circles) and the PPQ (squares) interactions, in the FTHFB 
approximation (empty symbols) and the averages with shape fluctuations (filled symbols). 
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FIG. 5: Results for ^ '^Dy with the Gogny force versus the temperature. Average values are presented as solid lines and filled 
circles. Selfconsistent FTHFB results as dashed lines and open circles. Upper row, panel (a), pairing energies for protons n and 
neutrons v. Panel (b), deformation parameters. Lower row, panel (c), excitation energies. Panel (d) specific heat calculated as 
defined in the text. Square symbols correspond to results using the Gogny force with pairing set to zero. In the temperature 
range from 1 to 1.5 MeV squares and circles are superimposed since above 1 MeV the solutions with and without pairing are 
the same. 
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FIG. 6: Same than Fig. I^for the nucleus ^^^Hg. In panel (d) the superposition of results with and without pairing starts at 
T=0.9 MeV. 
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FIG. 7: Total level densities in MeV ^ versus the excitation energy for FTHFB (empty symbols) and averages with shape 
fluctuations (solid symbols). 
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FIG. 8: Root mean squared radii, versus the temperature for FTHFB (empty symbols) and averages with shape fluctuations 
(soUd symbols). 



